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In this work, we explore the flavour changing decays Hi —>■ 6s in a general supersymmetric 
scenario. In these models, the flavour changing decays arise at loop-level but, originating 
from a dimension-four operator, do not decouple and may provide a first sign of new physics 
for heavy masses beyond collider reach. In the framework of the minimal supersymmetric 
extension of the Standard Model (MSSM), we find that the largest branching ratio of the 
lightest Higgs {Hi) is 0(10“®) after imposing present experimental constraints, while heavy 
Higgs states may still present branching ratios 0(10“®). In a more general supersymmetric 
scenario, where additional Higgs states may modify the Higgs mixings, the branching ratio 
BR(i7i —)■ 6s) can reach values 0(10“"^) , while heavy Higgses still remain at 0(10“®). 
Although these values are clearly out of reach for the LHC, a full study in a linear collider 
environment could be worth. 
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I. INTRODUCTION 

Since the discovery of ci scctlnr boson with insiss ■^126 GeV, nt the LHC in 2012 mm, the 
Standard Model (SM) picture may have been completed. Indeed, if this scalar particle corresponds 
to the SM Higgs boson, the SM could be the correct description of nature up to scales close to 
the Planck mass. So far, all the experimental evidence seems to be pointing in the direction of 
confirming that it is really the missing piece of the SM puzzle. Nevertheless the exploration of 
features of this particle is just beginning and further studies are needed to confirm its identity. 

From now on, our efforts to probe the SM and to search for physics beyond it may follow two 
complementary paths: i) push the energy frontier in the search for new particles and interactions 
and ii) increase the precision on the couplings of the first (so far) fundamental scalar ever discovered. 
Indirect searches, which would include the latter path, searches of rare processes or higher order 
corrections to low energy couplings, have been very successful in the past and have led to the 
discovery of new particles such as the third generation quarks. They have also been instrumental 
in exploring the scale of new physics beyond the collider reach. 

In the case of the Higgs boson, the study of its couplings can be the way to go. Scrutinizing 
non-standard Higgs couplings is a way to test the presence of additional scalar bosons even when 
their direct production is closed. In the models beyond the SM in which there exists more than 
one Higgs doublet, as in a two Higgs doublet model (2HDM)[3] or the MSSM [UlS], the couplings 
of the candidate for the discovered scalar boson may not be flavour diagonal and thus it can have 
flavour changing decays [6H231. Furthermore, these flavour changing couplings are dimensionless 
and therefore the effects of additional heavy particles may not decouple even when the masses of 
such particles are taken to infinity, providing a unique opportunity to find an assuring indirect 
evidence for such a high scale non-easily accessible in the near future. 

Even more, as we have no clue of the new energy scale (if any) associated with the new particles 
and interactions we are looking for, considering rare Higgs decays is a wise way to go. Thus, our 
proposal consists of searching for the Flavour-Changing (FC) Higgs decays. In the SM these FC 
decays do not exist and therefore their presence would undoubtedly signal the presence of new 
physics beyond the SM. 

As mentioned before, in the past, rare decays, including processes like 6 —)• sy and Bg —)• /x“, 
have been extensively used to search for new physics. Their precise experimental measurements 
were nseful for the exploration of the parameter space of different SM extensions. Likewise, FC 
Higgs decays are very useful to search for new physics since they are present in almost any extension 
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of the SM containing additional scalars that mix with the lighter Higgs. For instance they are 
unavoidable in type-II 2HDMs, like the MSSM, pseudo-dilaton models IMIES], models with extra 
dimensions 123 EH or composite Higgs models [28| . To be specihc, in this work we will explore 
a generic supersymmetric scenario, as well as the MSSM framework, both in the presence of non- 
minimal flavour structures. 

Our analysis will be focused on the process H ^ bs, since one would expect on general grounds 
that, among FC Higgs decays, those involving third generation particles whose Yukawa couplings 
are larger are the most experimentally accessible. Besides, loop-induced FC processes in the quark 
sector are typically larger by a factor (where and 02 are the coupling constants associated 

to the groups SU{3) and SU{2) respectively) for similar flavour changing entries in the lepton 
sector. 

The main goal in this analysis is to hnd out the largest FC branching ratios for the different 
Higgs states attainable in a general supersymmetric scenario. As we will show, in the MSSM, the 
decoupling of the heavy Higgses, enforced by present constraints, makes the FC branching ratio 
(BR) of the observed light-Higgs to be below the level of 10“®, while heavy Higgs states could 
reach ~ 10“^. In a more general supersymmetric standard model, the light Higgs BR could still 
reach values of ~ 10“^. Thus, these rare decays are very challenging for the LHC but they can be 
searched for at lepton colliders as the International Linear Collider (ILC). As will be shown, the 
branching ratio can reach a level of 10“^ under special circumstances in the generic SUSY case. 
These branching ratios could be reached at a linear collider. In fact, we can see that already at 
LEP, a limit of BR(Z —)• bs) < 0(10“^) was obtained with only 0(10®) Z-bosons in Ref. [29] . 
We can expect the larger statistics and improved experimental techniques to improve these limits 
(see also m)- In the case of the Higgs decays, we can produce between 0(10®) and 0(10®) Higgs 
bosons for niH = 125 GeV m, and therefore we can expect similar values for the lightest Higgs 
branching ratio. Thus, FC Higgs decays may provide an indirect hint for the existence of new 
physics at higher energies even when these higher scales are beyond the LHC or ILC reach. 

This paper is organized as follows: in section [HJ we introduce the framework in which the 
analysis will be carried out. Some compelling variations of it will also be addressed because of 
their interest when trying to observe FC processes. Still in this section, the theoretical expressions 


for the FC Higgs decays into down-type quarks will be provided. Section HI summarizes the 
latest experimental data corresponding to collider probes and indirect bounds from low-energy 
experiments. Finally, section [IV| collects the main features of the numerical analysis, concluding in 
section El with the main results. 
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II. HIGGS FLAVOUR CHANGING IN THE MSSM 


Our analysis is performed within a generic CP-violating MSSM framework and its extensions, 
in which the minimal Higgs sector is a type-II 2HDM, i.e. one of the two scalar doublets couples 
only to the up-type quarks at tree-level while the other couples only to down-type quarks. When 
electroweak symmetry breaking (EWSB) occurs, the neutral components of the two Higgs fields 
acquire vacuum expectation values (VEVs) and five physical Higgs states appear: three with neutral 
electric charge and two charged bosons. The two Higgs doublets can be parameterized as 




$2 






( 1 ) 


where 4 >? = {vi + 4>i + i Ui) , with vi = v cos 13, V 2 = v sin 13 and v ~ 246 GeV. At tree level the 
mass eigenstates (Hi) are CP eigenstates, but this situation changes once loop corrections are taken 
into account [32l - [38] . In the MSSM, the possibility of having CP-violating phases increases due to 
the growing number of complex parameters in the so-called soft SUSY breaking terms, and indeed 
these CP phases contribute at loop-level to Higgs masses and mixings. Consequently, weak-state 
fields and a) give rise to CP-mixed mass eigenstates (Hi) and these states are related through 
a unitary transformation represented by the 3x3 orthogonal mixing matrix O: 


4 >i — Hi, </>2 — 02i Hi, a — Ozi Hi. 


( 2 ) 


The mass eigenvalues will be obtained by means of diagonalising the mass squared matrix: 

O'^ ■ Mff ■ O = (3) 

To study Higgs flavour changing decays, it is helpful to introduce a convenient parametrization of 
the Higgs mixings. During the analysis, we will use 

S = f ) = ^31 / .s 

^ V cos 13 sin /? y sin j3 cos j3 ’ 

where di quantifies the distance of the lightest Higgs mixings from cos (3 and sin f3 and rji is directly 
related to the pseudoscalar content of Hi. 

In our analysis below, we will distinguish two different situations in regard to the Higgs mixing: 

• Full MSSM framework: Here, we consider the usual MSSM Higgs potential [5] which breaks 
the electroweak symmetry radiatively. The minimization of this potential gives us the Higgs 
masses and mixings. 
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Using of, + al, + al, = i, one may express the mixing angles Oai in terms of and ryi as 
follows 


Oil = cos /3 

\Jl — {51 + gl) cos^ /3 sin^ (3 + 5i sin^ (3 

£>21 = sin/3 

ijl — {5i+ gl) cos^ j3 sin^ j3 — 5i cos^ /3 

CI 31 = Vi cos /3 sin (3 . 


( 5 ) 


Then, the coupling of the lightest Higgs to a pair of massive vector bosons is given by 

ffnivv = cos /3 On + sin /3 £>21 = - ((5f + r/f) cos^ /3 sin^ (3 . (6) 

The current LHC Higgs data constrain g^^vv to be close to its SM value, gn^vv = 1- In fact, 
the present best-fit values and uncertainties are Ky = g^^yv = 1.15±0.08 if we assume there 
is no change in the Higgs total width, i.e. Kjj = T/z/TI^ = 1, and Kyy = Ky ■ Ky/Kn = 
1.28]’rQj5 if we allow for a change in the total decay width |39j . At present, the errors are still 
large, but requiring, for example, gniVV ^ 0.9, one needs to have (Sf+gi) cos^ /? sin^ (3 < 0.2. 
As we will show later, HR {Hi -^bs + sb) is directly proportional to the quantity ((5f -|- gf) 
which can be larger for larger tan/3 values while satisfying this constraint. On the other 
hand, large tan /3 values are constrained by the AB = 1 and AB = 2 processes such as 
b — )• S7, —)• 1 — B^^ mixing, etc. Taking into account all these constraint, we find 

that, HR {Hi -^bs + sb) can be as large as 10“® in an MSSM framework. 

• Generic supersymmetric SM: Given that no signs of supersymmetry have been found so far in 
collider experiments, and taking into account the strong constraints on the parameter space 
of minimal models, it is interesting to consider more general models. In fact, the situation 
could be different if we consider SUSY models beyond the MSSM which contain additional 
Higgs states. In this case, the Higgs mass eigenstates Hi are given by 

ai (t>a + O^i a -|- I] Ogi (7) 

0 = 1,2 / 3>4 

where ipp represent the additional CP-even and CP-odd Higgs states which can be charged 
or neutral under SU(2)x,. We note that only the SU(2)2,-charged CP-even states contribute 
to the tree-level g^.yy couplings and, due to the additional states, we generically have 
Oil + ^21 + ^31 < I- As in the MSSM framework, these couplings are constrained by the 
experimental results on Higgs decays, but in the presence of other Higgs states close to the 
125-GeV state, the mixing pattern could be different from that in the MSSM and and/or 
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r]i can be sizeable. In this case, one may treat di and rji as free parameters effectively. We 
find that, in this case, BR [Hi -^bs + sb) can be as large as 

Processes mediated by flavor changing neutral currents (FCNCs) involving down-type quarks 
have been largely studied in the context of 2HDM where significant contributions can be accom¬ 
modated due to the tan /3-enhancement of their Yukawa couplings. This type of processes are very 
useful for investigating the dynamics of quark-flavour mixing, especially the possible non-standard 
phenomena. Here, our main purpose is studying transitions such as Hi —>• bs, keeping always 
under control other processes that will impose additional experimental constraints, for instance, 
the B-meson decay Bg —)• and the mass difference AMb^ 


A. FC couplings 


It is well-known that, in the MSSM, the superpotential holomorphicity prevents the appearance 
of Higgs-boson FCNCs by coupling the Higgs-doublet superfield to the down-quark sector and 
<^2 to the up-quark sector. However, this property is violated when considering finite radiative 
(threshold) corrections due to soft SUSY-breaking interactions [i2lIM] . As a consequence, Higgs- 
mediated FCNCs show up at the one-loop level. The general effective Yukawa Lagrangian for 
down-type quarks may be simply written as m- 


—Cy — hd 


(l + A 


Vi + (/>! 

V2 






4,2 V2 + 4>2 




d^-|-H.c. 

( 8 ) 


where hd is the tree-level Yukawa matrix and dPi b to the weak eigenstates. After EWSB, 
this Lagrangian gives rise to the d-quarks mass terms and also to the Higgs-mediated FC terms. 
For the former, we have: 


'^mass 


Vl 

V2 


where 

basis: 


djj hd (1 -I- A^^^ hd A^^ +H-c. = ^ hd (1 + A^) d°+h.c. (9) 

+ (^^ 2 /^ 1 ) A^^ contains the loop corrections. Transforming the states to the mass 


4 = UUr-. dl =UidL = U'^YdL] (10) 

= ^rUR'i vPl = Ul , ( 11 ) 


^ In the presence of CP phases, limits from electric dipole moments (EDMs) should also be considered. However, 
in our scenario we have decoupled sfermions, and heavy Higgs masses are above the TeV. In these conditions, the 
main contributions to EDMs are due to scalar-pseudoscalar Higgs mixing from the two-loop Barr-Zee diagrams 
with Hi. As shown in snmi]. the relevant constraint on our couplings would be, |» 7 i| < (0.1 tan/3)^ which does 
not playing any relevant role in most of the considered parameter space. 
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we have: 

- ^mass = dR Mrf d^+b.C. with ^ {Kr} ^ hd [1 + Arf] V (12) 

where = diag (m^, m^) is the physical diagonal mass matrix for the down-type quarks. 

Using the flavour basis where = U'r = Ur = 1 and introducing = 1 -|- A^, we can relate the 
physical masses to the Yukawa couplings through the following expression: 

/o 

hd=—MrfVtR-i. (13) 

Vl ^ 

Using this expression in Eq. ([^, we obtain the FC effective Lagrangian for the interactions of the 
physical neutral Higgses with the down-type quarks [55l ES] : 


^FC 


g 

2 My/ 


Hi d (Md Pl + Md Pr) d 


(14) 


A simplified expression for the Higgs couplings can be obtained working in the single-Higgs-insertion 
(SHI) approximation |55j where: 


; 




(15) 


Ad = FO + = FO + tan/3 GO 


(16) 


For large tan/3 values, the term can be neglected and therefore R^ may be approximated as: 


Rd = 1-1- tan /3 G^ . 


(17) 


Then the Higgs couplings in the Lagrangian of Eq. (14) will be simplified to: 

Ou 


^Hidd 


^-VtR^iV + -^VtR^iGOV + iOsitanP 
cos P “ cos /3 a a 


— V'l'RT^V_ — 

.sin^/3 tan^/?. 


(18) 


^Hidd — {^H.dd) ■ (19) 

By noting = (l — ViR^^V) /tan/3, we observe that the size of flavour violation is 

dictated by the off-diagonal components of the matrix ViR^^V = (l-I-tan/3 Gj)) Therefore 
the Higgs couplings to the down-type quarks will be determined once G^ is known. 

The detailed expression for this quantity G° in terms of soft SUSY-breaking parameters can be 
found in Appendix [A} Observe that in this formalism, flavour violation from off-diagonal compo¬ 
nents of the sfermion mass matrices as well as its diagonal parts has been included. 
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B. FC Higgs decays 


For the computation of the decay width T[Hi bs + sb), we consider the relevant terms 


involving b and s quarks in Eq. (14). Introducing the effective couplings and 


VRi = we may write: 


^Hibs = -Hib ULi Pl + Ur, Pr 


s + H.c. 


This expression can be rewritten as: 

^Hibs — Hi b 


9i + wf 75 


s + H.c. 


( 20 ) 


( 21 ) 


with gf = (uLi + yRi)/‘2^ and gf = i{yLi - yRi)/‘^- Then, using Eq. (28) of [57|, the decay width 
can be obtained as: 

Nc ruHi 


r [Hi bs + sb) = 2 X 


lOvr 




Xb,Xs) kqcd X 


(1-xb- Xs) {\yLi? + \yR,?) - ^^/xbXs ?fie {yL,y*R^ 


( 22 ) 


where kqcd = 1 + 5. 67ay /vr including the QCD correction, Xf = rnj/rnj^, and A(l,a,5) 

(1 — a — b)'^ — Aab. Note that in yL.,_R. the masses involved are rrib^s = mb^s ("luj. 


III. EXPERIMENTAL CONSTRAINTS 

In our analyses, we use two main sets of experimental data: collider constraints on Higgs and 
supersymmetric particles and constraints on flavour changing processes. Collider constraints come 
mainly from CMS and ATLAS, the two general purpose experiments in the LHC that claimed the 
observation of a new 125 GeV particle in 2012 mm- Regarding indirect processes, we use the 
current flavour experimental data associated with B-meson decays and mass differences. 


A. Collider Constraints 

ATLAS and CMS are the two general purpose LHC experiments which provide the most accurate 
data concerning the Higgs boson and SUSY. In particular, we will consider the Higgs 77 signal, 
the rr-channel limits, and direct limits on supersymmetric particle masses. All these constraints 
have been already used in previous works |58( [59], so we refer to them for details. Here, we will 
summarize the basic requirements taken into account during the analysis. First, according to the 
experimental data so far, we require a diphoton signal in the range: 
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0.75 < < 1.55, (23) 

where fix is the signal strength for a Higgs decaying to X: fix = [cr{pp —)• H) x BR(H —)• 
X)]/[a{pp ^ H)sm X BR(H ^ X)sm]. 

On the other hand, we apply the limits set by CMS [60] and ATLAS m in the H —)• rr-channel. 
Specifically, we use the 95 % Confidence Level (CL) limits on the gluon-fusion and 6-associated 
Higgs boson production cross sections times the branching ratio into r pairs, presented in Fig. 4 in 
Ref. (GOj and Fig. 11 in Ref. ^Tj. In these analysis, extended searches for extra Higgs states have 
been carried out for masses up to 1 TeV at 95% CL. In our case, these limits will be imposed to 
all three neutral Higgs states: Hi, H 2 and H^. 

Finally, we must take into account direct bounds on SUSY masses. Taking into account that 
the effects of SUSY particles on Higgs couplings are non-decoupling, we can apply conservative 
limits on the masses. For the gluino, we set the mass limit at rrig ^1.4 TeV when the neutralino 
mass is below ~ 700 GeV, in agreement with the exclusion limits from ATLAS [62| and CMS m- 
The mass limits for the third generation squarks are taken from ATLAS data |62| as: ^ 650 

GeV when the neutralino mass is below 250 GeV, or — m^o ^ 175 GeV when it is nearly 
degenerate with the LSP. Finally, according to ATLAS searches |62) . the chargin mass limits are: 
m~± ^ 700 GeV for dominant decays into charged leptons, or m ± ^ 450 GeV when the decays 

Xl A1 

into weak bosons prevail. 


B. FC Constraints 

Apart from these data coming from the collider experiments, there is another kind of processes 
that can play a significant part in the search for SM extensions. As presented in the previous works 
[58l[59|, flavour constraints may be a powerful weapon to restrict the parameter space, especially 
the parameter tan/3, even in the absence of complex flavour structures beyond the SM Yukawa 
couplings. Therefore, for our analysis, we will make use of indirect bounds coming from B-meson 
decays and mass differences. In particular, we will consider: —)• fx^fi~, AMb^ and B —)■ 

In the case of the rare decay B^ —)• pHfji~, the latest experimental value for its branching fraction 
is the combined analysis of GMS and LHGb data at 7 TeV and 8 TeV, with integrated luminosities 
of 25 fb“^ and 3 fb“^ respectively [64] : 

BR(H0^/rV) = (2.81°;^) • 10-9 


(24) 
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Hence, our analysis bound at 2cr will be: 

BR ^ < 4.2 • 10“® (25) 

Another valuable flavour decay is R —)■ which becomes the most restrictive constraint for 

medium and low tan/3 values. Combining the BaBar, Belle and CLEO analysis, the world average 
value given by HFAG [65] is: 

BR (R ^ Xs-i) = (3.43 ± 0.21 ± 0.07) • 10"^ (26) 

For the opposite tan /3 regime, that is large tan /3 values, the main experimental result turns out 
to be the Bg-meson mass difference AMb^- The present experimental value is |65| : 

AMb, = (17.757 ± 0.021) ps"^ (27) 

We will require in our analysis: 

15.94ps"^ < AMb, < 19.83 ps"^ (28) 

where we included the theoretical error on /bsV^^ = 262 ± 10 |66) . 


IV. NUMERICAL ANALYSIS 


Before we present the results of our numerical analysis, we present first some approximate 
analytic expressions for Gjj, ^ and, finally, which are most relevant to the FC 

Higgs decay into b and s quarks. Note that the FC structure is common to the two situations 
under consideration: the full MSSM framework and the generic Supersymmetric SM. 

From Eq. (18) and using VtR-^G^V = (l - VtR-^v) / tan /3, we observe that the size of 
flavour violation is dictated by the off-diagonal components of the matrix V^R^^V which can be 
determined once G^ is known. 

As shown in Appendix [^ one may need the explicit forms of the flavour violating matrices 
(3 Mq jj and 6a^ to derive Gj). Assuming universality for the first two generations, we introduce 
the following flavour parametrization ^ 


M^ = 


^ p 0 0 ^ 

0 p Sx 
\0 6x 1 


M 


X3 ; 


= 


^ p 0 0 

0 P Sa, 
\0 1 




(29) 


We are assuming, for simplicity, symmetric Yukawa and trilinear matrices at tree-level. In general, (h„^a„)^^ and 


(h„^a„)^^ can be different from each other and complex. 
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where X = Q,U, D and then <5 Mq jj and (5a„ are given by 


^^Q,U,D — ^Q,U,D ~ ^Q,U,d'^ ^ 


SSyj, - 3.1, hijXi, 


(30) 


with = iTV = 5 (2P + 1) ^^^ 3,03 = S^r (h„^a„) = |(2p + 1) In 

the basis where the up-type Yukawa quarks are diagonal h„ = diag (?/„, Uc, yt), we have, 


5m\ = 


/ Pzd: 0 

0 ^ 
^ 3 


0 

<5x 


\ 


^ ^Vu 0 


M 


X 3 > 


(53.7/, — 


V 0 6 x -lip-l) J 


\ 


0 


^Vc 


^Au Vc 


V 


( 31 ) 


-f (d- 1) yt 


Inserting the above expression for (^Mg jj and 53u into Eqs. (A2) and (A3), we obtain 


/ 


g5~ 


e 0 0 

0 e (5623 


\ 


(32) 


0 (5632 6 + ?7 y 


where 6 , rj, 5 ez 2 and (5623 are parameters containing the main loop contributions of Eq. (A2) and 


Eq. (A3). Here, it is worth mentioning that EW corrections in those two equations have been 


neglected, while the down-type Yukawa couplings have been approximated as ~ ^M^Vt. The 
explicit forms of the diagonal entries e and rj are given in Appendix while the off-diagonal 
elements, which are the key ones for us as will be seen later, are given by the following expressions: 

|2 


6e23 = 5l 


+ 


y* M* Ml^K (mI,MI,\M3\^) + y* A*, 


\yt 


16 TT- 


;y*All(Ml,Mfj,\yi 


Sr 


Im 

16 

2 q!s Vij 3 Uh 
StT ^ 2*2 2 /s 


(33) 






)o ^ 


^632 = Sl 
+ Sr 


‘2as 

37r 


y*M*M^^K(^M^,Ml,\M3\ 

2Q:s f V22 Vs ^23 Vb 


<5a„ 


\yc 


IGtt^ 


y*Atl[M^,MlrM 


37 r y Y33 2/6 E22 Y33 2/s 

r 9rv 

- (p-1) 


y*M;Ml^K(^Ml,M^,\M3\- 


^^ p* M; M|,3 K IM 3 I" 


33 


(34) 


where 5q = Sy and Sjj = Sd ^ Sr. We note that there are four types of flavour-violating terms 
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proportional to 6 l, Sa^-, Sr, and {p — 1) Then, we obtain 


R-l ~ 


( (l+£tan/3)(l+(e+r)) tan/3)—(5e23<?<;32 tan^ fi 

Det(Rd) 

0 
0 


0 


0 


from 


/ 


Rd = 1 + tan /3 


(l+e tan/3)(l+(e+T;) tan/3) (5e23(l+e tan /3) tan /3 

Det(Rtj) Det(Rd) 

'5e32(l+etan/3) tan/3 (l+etan/3)^ 

Det(R(j) Det(Rtj) 


\ 


(35) 


1 + e tan /3 0 0 

0 1 + e tan (3 6623 tan (3 

y 0 (5e32tan/3 1 + (e + r/) tan/3 y 

As we will see, the most relevant flavour-violating matrix element in the FC Higgs decay Hi 


(36) 


bs 


is (vtR^V) . In principle, this matrix element (V^R^^Vj is: 
\ / 32.23 ' ' 32 


*J =1 


E [V?3 (itj')ii ’^“1 + '"a K')2S ‘)32 


2=1 


Using Eq. (35) for R^ and taking into account the unitarity of the CKM matrix, we can write: 


^ [v:3 V,2] = (V^3 Vi2 + V^gV: 


2=1 


(1 + etan,d) r/tan/3 (3e235e32 tan^/3 

’ 22 j- , /-r. \ -V 33 V12- 


Det (Rd) 


Det (Rd) 


(38) 

In this expression, we can neglect all the terms proportional to VJ 3 V 12 ~ 8 x 10“^ with respect 
to V 23 V 22 ~ 4 X 10“^, even for the last term proportional to 5e23Se32, which, as can be seen from 
Eqs. ( [33|34l ) and ( |A6|A7D , is, for sizeable mass insertion (MI), of the same order as e x p. 


Then, in Eq. (37), if we have similar values of the off-diagonal elements and 

the former can also be neglected with respect to later, being suppressed by an additional | V 23 V 32 I ~ 
2.10“^. Therefore, in the presence of sizeable mass insertions 6 l,r > V 23 V 22 , we have 6632 > 
p X V 23 V 22 and then we can safely take, 

(I -|- etan^d) p tan/3 


V^R-^V 


32 


V^3V22 


(IIS ')32 + 

Repeating the same exercise with we obtain, 


Det (Rd) 




V^2V33 


(I -|- etan/3) p tan/3 
Det (Rd) 


(R.-‘)32 ■ P") 


(R.-‘)23 ■ 


3 


Please note that our definition of 6 au in Eq (311 makes it different from the Slr usually defined in the literature. 
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The study of these matrix elements is very interesting because of their dependence on 5 l, Sr and 


6 au- Looking at Eq. (34) for 6€32 and comparing the 6 l and Sr contributions, we can see that the 


Sr term is suppressed by the difference (jjs/Vb — ^ 2 ^yb/y^ — —0.013 with respect to the Sr term. 
Therefore, these matrix elements, and especially will have very different values for these 

two mass insertions. Also from these equations, it is evident that, if Sr — Sa^j we would obtain 
similar results from both types of mass insertions to S€ 23 , but its contributions to Se 32 would be 
suppressed by an additional factor {yc/yt)‘^- In any case, the Sa^ contributions to S €23 and Se 32 are 
always smaller than the Sr contributions. For this reason, we will only consider the cases of the 
Sr and Sr insertions 

All the results of the numerical analysis presented in the following sections and the corresponding 
figures are done with the CPsuperH2.3 code [5711671168] . 

A. Full MSSM framework 

In the full MSSM framework, the effective FC Higgs couplings are given by: 


yLi 


- ^bL 
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^Hibs ~ 


rUg R 

yRi = —^Hy,s 


m.g 


and the decay width, Eq. (22), is: 


T {Hi ^bs + sb) 
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Sttu^ 
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li 


02^ 


VtR-^V 
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COS fi sin /3 

'■(vtRj'V 


+ tan/3 — 


03i 
s\n^ fS 


2 m- 

4- —^ 
' 2 

mt 


23 


In the case of Sr insertions, from the discussion in the previous section, we observe <C Ii/lJ 

due to the ms/rrih suppression with . In the presence of Sr insertions 

\ / 32 N / 23 

the situation is more involved and both terms must be considered. 

Now, considering the total decay widths of the Higgs bosons, we will obtain the corresponding 
branching ratios. In the case of the lightest Higgs, its total decay width is dominated by the decay 
into two 6 -quarks, two VF-bosons and two r-leptons. Thus: 

‘ + 3 ) ^ (^'1 + ^3^) + Ips^ (021 + (44) 


Thi = 


mH^rnl 

Svrn^ 
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where in the large-tan {3 limit we have used {Oh/ cos^ f3 + tan^ {3 O^) — tan^ [3 {Oh + Oh) and 

(sin/3 O 21 -|-cos/?Oii)^ ~ (O 21 + dii/tan/3)^. Ips in the second term refers to the phase-space 
integral in the Higgs decay into two W bosons m and can be approximated by Ips — 6.7 x 10 
when niHi ~ 1^5 GeV. Then, the branching ratio will be: 


BR {Hi -^bs + sb) = (V^R-^V) 


32 


3 KQCD 


COS 0 sin 


O2: 




031 

sin (3 cos (3 


(3 nQCD + tan2 f3 {Oh + Oh) + Ips^ (o,i + 


2 • 


(45) 


For the heavier Higgses, and for tan (3 ^ 30, the total decay width is dominated by the bottom 
and tan widths: 






'f^QCD + —§ I tan^ j3 foh + oh 

mb J 


(46) 


using {oh/ cos^ {3 + tan^ /3 O 32 ) tan^ /3 {Oh + Oh) and the branching ratio is: 


BR {H 2 ^bs + sb) =3 kqcd I (vtR-iy) 


32 


( 012 _ 022 y 

y cos ^ sin (3 ) 




/ 032 ^ ' 

Y sin /3 cos /3 J 


(47) 


3 Kqcd + ^ j tan^ (3 {Oh + ^ 32 ) 

In a previous work [59], we showed that the latest LHC data for the Higgs signal in the diphoton 
channel strongly constrains the Higgs mixing within a general CP-violating model. In particular, 
if we consider the lightest Higgs as the recently discovered boson at 126 GeV, its mixing conditions 
are{Oh + Oh) 1 / tan^ /3 and O 21 1. Additionally, using the parametrization presented in 
Eq. @) , we have: 


BR {Hi ^ 6 s + sb) = 


V 


Ri'v) 


32 


3 kqcd (<^i +??l) 


3 Kqcd + ^ 


"If. 


(48) 


Still following |59j , it was also showed that the diphoton condition establishes for the heavier Higgs 
mixings that: (O^j -I- Ol/) ^ 1 and © 2 * ^1/ tan^, i = 2,3. Hence: 


BR(R2 ^^s + s6 ) ~ |(vfR^^V 


3 Kqcd 


(49) 


' *^9 O 1 

3 kqcd+ 

where we have considered tan/3 > 3 and, therefore, l/cos/3 ~ tan/3 and sin/3 ~ 1 in a good 


approximation. 


1. Left-handed (L) insertion 

First, we analyze the case with Sp 0 and 6p = 0: 

-5e32(L) ^ ~-3-10-36 c (50) 
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BR(N,->bs) 
10"'* r 



BR(Ni-»bs) 


laniP) 


10"^ 
10 "® 
10"® 
10"^ 
10-® 
10"® 


(<5iP 

FIG. 1: A full MSSM framework with LL insertion with 6l ^ 0 and 6a^ = = 0: The upper frames 
show the dependence of the estimated branching ratios for Hi 2 bs + sb on tan/3. The lower-left frame 
is for the dependence of B{Hi —>■ bs) on 6f and the lower-right frame for the Sf dependence on tan/3. Blue 
(dark) points satisfy all the constraints considered, while red (light) points violate one or several of these 
constraints. 



BR( H2 “* bs) 
10-= 


taii(p) 




where we can see that Se^ 2 {L) has a non-decoupling behaviour as it depends only on ratios of 
sparticle masses and we have used ~ 5 TeV, M3 ~ 7 TeV and ~ 6 TeV. Consequently, 

for the maximum values of tan /3 and 6l considered during the scan, 6l ~ 0.5 and tan j3 ~ 60, we 
would obtain: 


(rF) 


32(1,) 


<5e32(L) (1 + etan,d) tan/3 3-10“^ x 1.5 x 60 
Det (Rd) 4 


0.03 


(51) 


where e, Eq. (A6), and Det(Rd) for the masses specified above take the values e ~ 0.01 and 
Det(Rd) ~ 4.7 . Therefore: 

2 Skqcd X {6l + r]f) 


BR {Hi ^bs + s&)“ 
BR {H2 -^bs + sb) 
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(VtR-iV) 


(VtR-^V) 
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(skqcd + ^) + Ip S 
3 kqcd 
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o I mi 
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-3 


5 • 10-^ (df + ) (52) 


(53) 


Notice that, effectively, there is no bound from low-energy FC processes on this MI for such heavy gluinos and 
squarks. 
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Fig.[^shows the results of our scans. Blue points are those which satisfy the whole set of constraints 
while red points are excluded because of the violation of one or more of them. In the upper frames 
of Fig. we represent the branching ratios for Hi and H 2 versus tan /3 (for all considerations, 
will be equivalent to H 2 given that they are nearly degenerated for the range of masses considered 
here). As can be seen in these two plots, before applying the low-energy constraints, the branching 
ratio grows with tan/3 (red points). However, the final result is very different from this when we 
impose experimental limits from B-mesons (blue points). Whereas for heavy Higgses few points 
become excluded, in the case of the lightest Higgs the effect is more notable. Indeed, looking at the 
upper-left frame of Fig. [l] we can say that the tendency is completely opposite and the branching 
ratio decreases for tan/3 > 20. 

This behaviour is mainly due to the —?• constraint. This branching ratio is given 


in Eq. (BI3) with Cs, Eq. (BI5), and Cp, Eq. (BI6), containing the SUSY contributions. The 


dominant contributions come from the heavy Higgses and we have: 


tan2/3 
oc 2 —2— 
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Hi 


(vtRjW) 


32 


Oii- 
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2i 


tan^ (3 
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tan^ /3 
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H, 


(r.-‘) 


32 


(54) 


with j = 2,3. This dependence on tam/3 (the matrix element (R//^)^^ carries an additional tan/3) 
and the heavy Higgs masses explains why this decay provides such a restrictive constraint for 
relatively small and medium-to-large values of tan /3. 

In fact, BR(Ffi —)• hs) is suppressed at medium and large tan/3 values while BR{H 2 —)• bs) is 
not. This is due to fact that BR(Ffi —)> bs) is proportional to (3i and r/i, which in the MSSM are of 
order Then, B-meson constraints are more restrictive for light mH 2 which also correspond 

to the largest BR(Ffi —?• bs). On the other hand BR(Ff2 bs) is independent of (5i or rji and the 
H 2 mass. The B-meson constraints are not effective here if mH 2 is large enough and therefore we 
can reach large branching ratios for large tan/3 values. Eurthermore, this branching ratio saturates 
for medium-to-large values of tan /3 when both the EC decay width and the total decay width have 
the same tan /3 dependence. 

Moreover, the lower-left frame of Eig. shows the branching ratio for Hi versus . As seen in 


Eq. (52), for Sf ^ 0.02 the branching ratio is of the order I0“®. However, the implementation of 
B-meson constraints here reduces this value by more than one order of magnitude. 

Einally, in the lower-right frame of Eig. we present the 5i dependence on tan /3. We observe 
that, before imposing B-meson constraints, larger values of (5^) are possible when tan/3 grows. 
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BR(H,->bs) 


- tan(P) 


BR(/y2->bs) 



„tan(jS) 


FIG. 2: A full MSSM framework with RR insertion with Sji ^ 0 and = 0: The left frame shows 

the dependence of BR(i7i — bs) on tan/3 and the right frame shows the dependence of BR(i72 —5^ bs) on 
tan p. Again, blue (dark) points satisfy all the constraints considered, while red (light) points violate one 
or several of these constraints. 

However, B-meson constraints become very effective for tan /3 > 17 as shown by the blue points 
in this plot. This could be understood by noting that the Higgs contributions to the AB = 1 
and AB = 2 processes are inversely proportional to the heavy Higgs mass squared. Therefore, to 
suppress these processes for large tan /3, large Higgs masses are required. Accordingly, we expect 
smaller |(5i| as tan/3 grows since, as we have seen, |(5i| oc in the MSSM. 


2. Right-handed (R) insertion 


Now we consider the case, and 6l = = 0) we have from Eqs. (34) and (50): 

'^hVs ^h^Vb \ <Je32(L) 




~ 0.013 X 3 X 10"^ ~ 5 • 10"^ 5r , (55) 


-5 






(56) 


In this case, the value of the off-diagonal element (R^ )32 is of the order of 10 ^ and this implies 


that the contributions V23V22 in Eq. (38) and (R^^)23 in Eq. (37) or yR. in Eq. (42), can be 
important. Thus: 
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where we used the same values for the parameters as in Eq. (51) and (5e23(R) is now enhanced by 
a factor yb/Vs'- 


Se23(R) ~ y* M; Ml^K (m^,MI,\M 3 \^) -O.lSn. 

OTT 


Therefore, we obtain, 


(vtR-vr“ 


-3 


32{R) 


Thus, for 5r = 0.5, the branching ratios are: 


max 

(fl) 


BR (i?! —>■ 6s + s6) 

BR {H2 ^bs + s6)“" 




4{6r + 1 ) • 10 


3 KQCO + ? 7 l ) 


(3«:qco + ^) +Ips^ 
^ 3 kqcd 


o I 

3kqcd + 


2 • 10 
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( 58 ) 


(59) 


1.5- 10-'^(<5f+ r?fX60) 


( 61 ) 


The results of our scans for this case are shown in Fig. As before, these results are in 
agreement with the numerical values if B-meson constraints are not taken into account. Once they 
are incorporated into the analysis, the lightest Higgs branching ratio is reduced by more than one 


order of magnitude. Also, taking into account, from Eq. (57), that both the MI and Ml-independent 


contributions are of the same order only for large 6r, the BR is completely independent of 5r. 


B. Generic supersymmetric SM 


After computing these branching ratios in the MSSM framework, we now perform our analysis 
in a generic supersymmetric model. Therefore, we present here a model-independent analysis, 
meaning that, in fact, we consider a generic Higgs mixing matrix with possible additional Higgs 
states. In this case, we have Oh + O 21 + O^i < 1 and the parameters and/or 771 entering 
T{Hi -^bs + sb) can be sizeable. 


The expressions for the decay widths and branching ratios, Eqs. (43-49), are still valid in the 


generic supersymmetric scenario. The main difference now is that the parameters (5i and r]i are 
only constrained by experimental results on Higgs decays and low-energy FCNC processes. Notice, 
however, that the flavour-changing entries in do not change in the two models. 

In Fig. we show the EC branching ratios of Hi and H 2 for <5^/0 and 5au = = 0 

in the generic supersymmetric SM. These figures can be compared with Fig. which shows the 
corresponding branching ratios in the MSSM framework. 
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BR(H,-»bs) BR(H2-.bs) 



FIG. 3: A generic supersymmetric SM with L insertion with Sl ^ 0 and <5^4^ = di? = 0: The upper frames 
show the dependence of the estimated branching ratios for Hi 2 —>■ bs + sb on tan/3. The lower-left frame 
is for the dependence of B{Hi —>■ bs) on and the lower-right frame for the Sf dependence on tan/3. Blue 
(dark) points satisfy all the constraints considered, while red (light) points violate one or several of these 
constraints. 


In the MSSM framework, the mixing angles are obtained through a minimization of the scalar 
potential and both <5i and rji are of the order of v '^from the diagonalization of the neutral Higgs 
mass matrix. In the generic supersymmetric scenario we treat 5i and r]i as free parameters that do 
not depend a priori on the ratio but are only constrained by the different experimentally 

measured Higgs branching ratios and B-meson constraints. This is why BR(ifi —?■ hs) in Fig[^is two 
orders of magnitude smaller that the largest possible value in Figure The different distribution of 
the points allowed by B-meson constraints is due to the same reason. In the MSSM scenario, Higgs 
flavour changing processes are mediated by the heavy Higgses and therefore are only important 
for light mH 2 which, as we have seen, correspond also to the largest (5i and t/i and therefore to the 
largest branching ratios. In the generic supersymmetric SM, it is, in principle, possible to have a 
large 5i with a heavy H 2 and therefore the B-meson constraints are not so efficient. 

On the other hand, the FC decays of the heavy Higgses ^^2,3 are independent of the values of 
and r]i as can be seen in Eq. (49). Therefore, the upper-right frames of Figs. and are very 
similar and we obtain very similar results for BR(Ff2 bs) in the MSSM framework and in the 
generic supersymmetric SM. 
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Also, as shown in the lower-right plot in Fig. the allowed values of (5i)^ are completely 
independent of tan/3 as the B-meson constraints, which depending on (tan/ 3 )"'/m |^2 can always 
be satisfied by adjusting conveniently the value of In this case, the upper limit for (di)^ 

is fixed by the Hi —)■ 77 decay which as shown in [59], requires (Oji + O^i) ~ l/tan^/3 and 
O 21 ~ 1 — 1/ tan^ /3. Using the definition of (5i and in the limit O 31 <C 1, with the above constraints, 
(di)^ < 0.17 for tan/3 > 10, as we see numerically in this plot. 

In summary, the main difference in the generic supersymmetric SM is that BR(F7i — bs) could 
reach a value of ~ 10“^ consistently with present experimental constraints. This value is still too 
small to be observed in the large background of a hadron collider, but it could be tested in a 
leptonic linear collider in the near future. 

V. CONCLUSIONS 

In this paper, we have analyzed the FC Higgs decay Hi —)• bs for the different Higgs states in 
both an MSSM scenario and a more general supersymmetric SM framework. The importance of 
this observable is that effects of heavy particles do not decouple and may provide a first sign of 
new physics for heavy supersymmetric masses beyond collider reach. Before we carried out our 
numerical analysis, we derived approximated analytic expressions for the off-diagonal entries of 
(ytR-iy) which dictate the size of flavour violation. In an MSSM framework we showed that, 
even in the presence of large off-diagonal flavour entries in the sfermion mass matrices, for the light 
Higgs, BR(i7i —)• bs) ^ 10“® consistently with present experimental constraints, while for heavy 
Higgs states BR(i 72,3 —?■ bs) can still be ~ 10“^. In a more general supersymmetric scenario, where 
we allowed for non-minimal Higgs mixings, the branching ratio BR(F7i —)• bs) can reach values 
~ 0(10“^), while BR(F 72,3 bs) remain of the order of ~ 10“^. We find that the results of the 
numerical analysis are well in accord with the estimations made using the approximated analytic 
expression for ^yfR/j’^y^ Although these small branching ratios are clearly out of reach for the 
LHC due to the very large 6 -quark background, a full study in a linear collider environment could 
still be worth pursuing. 
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Appendix A: FC Higgs couplings 


In this appendix, we present the explicit expression for associated with the FCNC Higgs 


couplings in Eq. (14), 
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where the loops functions are given by: 
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From here, the elements e, 6, used in the matrix in Eq. (32) are. 
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Appendix B: B-physics constraints 


The main FC processes associated with B-Mesons that we consider in our analysis are AMb^ 
and Bg —>■ , although other constraints like BR(i? —>■ are also included. 

In the case of AMb^, we use the expression in |55j . given by: 


AMb^ = 2 
where the SUSY contribution is [55] : 
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with double-penguin and box diagrams. 

The Wilson coefficients ^SRR{dp)^ and related to the SUSY 


contribution of the B^-meson mass difference in Eq. (B2) are 
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The couplings of the charged Higgses and Goldstone bosons to fermions, Eq. (3.38) of Ref. [56] . 
are given by: 
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and in the large tan/3 limit, we have: 
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Then main contributions: one associated with box diagrams for two Hj^ 

and another for box diagrams. From Eqs. (4.4) and (4.5) in Ref. [69], we have: 
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where DQ{a,b,c,d) and D 2 {a,b,c,d) are the corresponding loop functions which can be found in 
Ref. [69] 

The decay —)• is described by the effective Hamiltonian, 
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where the relevant operators are Os = (gPub) (pfJ-), Op = ^^rrib (qPpb) and Oiq = 

iqi^PLb) 

Neglecting the non-holomorphic vertices on the leptonic sector as well as the contributions 
proportional to the lighter quark masses ^, the branching ratio is given by 
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where Tp, is the total lifetime of the Bg meson and the form factors are: 
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where Cio is the leading SM contribution, and and — ~ ^3* 

Higgs couplings to the charged leptons. 
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